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Supernova neutrino halo and the suppression of self-induced flavor conversion 
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Neutrinos streaming from a supernova core occasionally scatter in the envelope, producing a small 
"neutrino halo" with a much broader angle distribution than the primary flux originating directly 
from the core. Cherry et al. have recently pointed out that, during the accretion phase, the halo 
actually dominates neutrino-neutrino refraction at distances exceeding some 100 km. However, 
the multiangle matter effect (which increases if the angle distribution is broader) still appears to 
suppress self-induced flavor conversion during the accretion phase. 

PACS numbers: 14.60.Pq, 97.60.Bw 



I. INTRODUCTION 

Neutrino-neutrino refraction is responsible for the in- 
triguing effect of self-induced flavor conversion (TMHl that 
can occur in the neutrino flux streaming from a supernova 
(SN) core In this context, the angular neutrino distri- 
bution plays a crucial role. The current-current structure 
of low-energy weak interactions implies that the interac- 
tion energy between two relativistic particles involves a 
factor (1 — cos 6*) where 9 is their relative direction of 
motion. 

If the neutrino-emitting region of a supernova core 
("neutrino sphere") has radius R, then at distances 
r ^> R a typical neutrino-neutrino angle is 9 ~ R/r and 
(1 — cos6>) oc (R/r) 2 . The geometric flux dilution pro- 
vides another factor (R/r) 2 , leading to an overall (R/r) 4 
decrease of the neutrino- neutrino interaction energy [6j . 

In a recent paper, Cherry et al. pointed out that this 
picture is not complete because neutrinos suffer residual 
collisions on their way out 7]. Every layer of matter 
above the neutrino sphere is a secondary source, produc- 
ing a wide-angle "halo" for the forward-peaked primary 
flux. While the halo flux is small, its broad angular dis- 
tribution allows it to dominate the neutrino-neutrino in- 
teraction energy. 

We illustrate the halo with a numerical example, the 
280 ms postbounce snapshot of a spherical 15 Mq model. 
We recently used it as our benchmark case to study mul- 
tiangle suppression of self-induced flavor conversion yj. 
In Fig. [T] we show the angular dependence of the inten- 
sity 1 for the v e radiation field, normalized to the forward 
direction, measured at the radial distances 300, 1000, 
3000, and 10,000 km. (The angular distributions become 
noisy for 9 > tt/2, where they are currently not well pro- 
vided by our simulations.) The core and halo fluxes are 
two distinct components, the latter so small that it is not 
visible on a linear plot. If we use 9 C ~ 0.1 as the edge 



of the core distribution for the 300 km case, we infer a 
radius of R ~ 30 km for the region where neutrinos begin 
to stream almost freely. At larger distances, the angular 
scales are squeezed by a factor R/r. 

The impact on neutrino-neutrino refraction is illus- 
trated by the weak potential felt by a radially moving 
neutrino. In terms of the zenith angle 9 of the intensity 
T(9) we need the quantity 



f£d9 sin (1 - cos 9)1(9) 
^ d9 sin 91(9) 



(1) 



For small 9, the integrand in the numerator expands as 
(1 — cos 9) sin 9 = 9 3 /2, allowing the halo to contribute 
significantly at larger distances where the distributions 
are more squeezed. At 1000 km, for example, (1 — cos#) 
is almost a factor of 10 larger than it would be without 
the halo, in agreement with Cherry et al. [7[. The halo 
parts of the functions in Fig. [T] decrease roughly as 9~ 3 
and we find this implies that (1 — cos 9) decreases roughly 



as r instead of : 



The neutrino-neutrino interaction 



energy then decreases roughly as r -3 instead of r~ 4 . 

The halo is important during the accretion phase when 
there is enough matter for the primary flux to scatter 
on 0- However, the same high matter density tends to 
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With "intenstiy" I we mean the quantity "neutrinos per unit 
area per unit time per unit solid angle," integrated over the en- 
ergy spectrum. 



FIG. 1: Intensity for the v e radiation field of our numerical 
model, normalized to the forward direction, measured at the 
radial distances 300, 1000, 3000, and 10,000 km (right to left). 
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suppress the self- induced flavor instability [9(. During 
the early accretion phase, self-induced flavor conversions 
were found to be typically suppressed [1, H3] • Does the 
halo change these conclusions? Its importance derives 
from its broad angle distribution, which also makes it 
susceptible to the multiangle matter suppression. 

To investigate this question, we study in Sec. |TT] the 
properties of the neutrino halo. In Sec. IIIII we perform a 
stability analysis along the lines of our recent study of an 
accretion-phase model 0, and we conclude in Sec. IIVI 



II. NEUTRINO HALO 
A. Numerical angle distribution 
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FIG. 3: Energy spectrum of core and halo components. His- 
togram: Numerical results. Smooth lines: Thermal spectra 
as described in the text with the same T. 



Neutrinos stream almost freely and therefore, at larger 
distances, the angular distributions are simply squeezed 
to smaller angular scales (Fig. [IJ . This behavior applies 
also to the halo flux which, once produced, streams al- 
most freely. At larger distances, the halo distribution 
primarily gains at the edges: the newly available angular 
modes get populated by scattering. 

The numerical angular distribution at the largest avail- 
able radius essentially holds all the required information. 
The decreasing wings of 1(9) look like power laws and 
we can fit the entire distribution by 



2fit(0) 



0.9994 



[1 + (6I/0.0029) 4 - 5 ] 2 
0.0006 



[l + ^/O.Ol) 1 - 



1/5 



(2) 



B. Energy distribution 

The v e flux spectrum emitted from the core roughly 
follows a thermal Maxwell-Boltzmann form f c (E) oc 
E 2 e~ E / T '. In Fig. [3] we show as a histogram (blue) the 
numerical spectrum of the core component (measured in 
the forward direction) together with a thermal fit with 
T = 4.8 MeV (average energy 14.4 MeV). 

The halo component arises from scattering on nuclei 
with a cross section proportional to E 2 . Therefore, the 
halo spectrum should be fh(E) oc E i e^ E / T with the 
same T . In Fig. [3] we show as a red histogram the nu- 
merical halo spectrum (measured at a very large angle) 
together with such a fit. In the spirit of an overall con- 
sistency check, we find excellent agreement. 



Neutrino-neutrino refraction 



We show this function overlaid with the 10 4 km case in 
Fig. [5J For the other species, the situation is analogous 
with slightly different parameters. 





10° 




10- 2 






to 


10- 4 


c: 




£ 




J: 


10~ 6 




10- 8 




io- 10 




10 -ii 10~ 1 
Zenith Angle 9 [rad] 



FIG. 2: The 10 4 km numerical case (solid line) overlaid with 
the analytic fit (dashed line) of Eq. ([2]). We also indicate 
the power-law behavior estimated from an analytic argument 
(Sec.|ITD|. 



The impact on neutrino-neutrino refraction is quanti- 
fied by (1 — cos 9) as defined in Eq. (JTJ) . Motivated by the 
numerical examples, we first consider truncated power- 
law intensity distributions of the form 



1(6) = 



1 for 9 < 9 C 

(9 c /9)p otherwise. 



(3) 



We ask for the asymptotic behavior of (1 — cos 9) at large 
distance, corresponding to 9 C <C 1. The integrand in 
Eq. (P) expands as sin6>(l - cosO) 9 3 /2. If p > 4, 
the integral is dominated by small angles and we find the 
result shown in Table |U In other words, if 1(9) falls off 
fast enough, we recover the classic r -2 scaling, where we 
have used that 9 C ~ R/r. 

For p < 4 we can no longer extend the upper integra- 
tion limit to oo and can no longer expand the integrands. 
With Mathematica we find analytic results with coeffi- 
cients for 9 C <C 1 that are given in Table Q] In particular, 
for p — 3 the scaling is linear in 9 C . 

Next we consider the analytic fit of Eq. ^ and show 
(1 — cos 9) in Fig. 2] We integrate up to 9 = it for each ra- 
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TABLE I: Average (1 - cos6>) for 9 C < 1 and different p. 



Power p 



(1 — cos ( 



p > 4 

4 

3 

2 

1 




4 c 

P-2 £ 
p — 4 4 

(1.1897 - 2 log 6 C 

0.54033 C 

1.5788 
1.9929 - 2 log 9 C 
0.61712 
1 



dius, which means that we smoothly interpolate between 
small and large distances. At small radii (inside of the 
SN core) the distribution is isotropic and (1 — cos6>) = 1. 
For r ;$ 200 km we find approximately the naive r~ 2 
scaling. At larger distances, the halo becomes important 
and the scaling turns approximately to r _1 . If the halo 
distribution would scale as 9~ 3 we would expect precisely 
r -1 , but our fit corresponds to 1(9) oc Q~ 2 m , explaining 
the less steep variation. 

In Fig. [5] we show the enhancement of (1 — cos 6) caused 
by the halo, where we attribute the first part of the ana- 
lytic fit of Eq. © to the core. The enhancement caused 
by the halo scales almost linearly at large distances. At 
r ~ 1000 km the enhancement is about a factor of 8, 
roughly in agreement with Cherry et al. [7j|. 



D. Analytic halo estimate 

For an analytic estimate, we consider a total neutrino 
production rate $ = L/(E) emitted from a pointlike 
source (neutrino luminosity L). It traverses a spheri- 
cal matter distribution which we model as a decreasing 
power law of the form n(r) = nu(R/f) m . We assume 
that multiple scatterings can be neglected. Every spheri- 
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FIG. 4: Average (1 — cos#) for our analytic angle distribution 
of Eq. ^ , representing the 280 ms numerical model. 
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FIG. 5: Same as Fig.[4l now showing the enhancement caused 
by the halo, i.e. the ratio total/core contribution, where the 
core flux is attributed to the first term in Eq. ([2]). 



cal shell is a neutrino source from scattering the primary 
flux. For the scattering cross section, we assume the form 
da/dft = (7scatt(l + acos#)/47r, where -d is the scattering 
angle. Elementary geometry (see the Appendix) reveals 
that at distance r we expect a scattering flux (halo) of 



l(r, 



(4tt) 2 R \r sin 

d&{\ + a costf) sin m i9. 



(4) 



For small angles, we find a power law I(r,Q) ot Q~( m+1 \ 
whereas the large-angle and backward flux depend on the 
detailed cross-section angular dependence. 

In our numerical example, the density outside the 
shock wave (at about 70 km) decreases roughly as r~ 135 
out to about 5000 km. With m = 1.35 we expect the 
halo to vary as 0~ 2 - 35 as indicated in Fig. [5] While it 
would not provide a good global fit, it looks excellent for 
about the first decade of angles of the halo. 

Based on the analytic expression for I(r, 9) we can es- 
timate the inward-going flux. It is very small compared 
to the outward-going core flux, but its contribution to 
refraction need not be small. We consider a region at 
sufficiently large radius where neutrino-neutrino refrac- 
tion is dominated by the halo and ask which fraction is 
caused by inward-going neutrinos, i.e. which fraction of 
the integral in Eq. (JXJ) is contributed by ir/2 < 9 < ir. For 
sufficiently small power-law index m this fraction does 
not depend on the radius r (to lowest order). Coherent 
neutrino-nucleus scattering is forward peaked and we use 
a = 1. For m = 1 the contribution of the backward flux 
is 25%, for m = 2 approximately 16%. 



III. STABILITY ANALYSIS 

All previous studies of collective flavor oscillations of 
SN neutrinos used the free-streaming approximation: it 
was assumed that collisions play no role in the relevant 
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stellar region. On the other hand, the halo dominates 
neutrino-neutrino refraction at larger radii, so in some 
sense the free-streaming assumption gets worse with dis- 
tance, not better. Therefore, in order to understand col- 
lective flavor conversion in this region, one may have to 
rethink the overall approach [7|- 

On the other hand, the multiangle matter effect may 
still suppress the onset of collective flavor conversions, at 
least in some accretion-phase models. A linearized sta- 
bility analysis would still be possible and self-consistent. 
If self-induced flavor conversions have not occurred up 
to some radius and if the matter effect is large in that 
region, then neutrinos continue to be in weak-interaction 
eigenstates up to small corrections caused by oscillations 
with the small in-matter effective mixing angle. There- 
fore, we can use the core+halo flux at that radius as an 
inner boundary condition for the subsequent evolution. 
In other words, collisions at smaller radii do not change 
the conceptual approach, they only provide a broader 
angular distribution than would have existed otherwise. 

A greater problem is how to deal with the inward-going 
flux which, for our conditions, contributes some 20% to 
neutrino-neutrino refraction. We will simply neglect it 
because it cannot be incorporated self-consistently in an 
approach based on outward-only streaming neutrinos. If 
the system is stable and stationary, then the picture re- 
mains self-consistent because neutrinos remain in flavor 
eigenstates, including the backward-moving ones. 

In this spirit we repeat our linearized flavor stability 
analysis for our 280 ms snapshot of a 15 M accretion- 
phase model Q in a simplified form. We model the an- 
gular distribution according to our analytic fit, Eq. ([2]l. 
which we use at all distances in the form 

X(r,0) oc 2fi t (0 r/10 4 km), (5) 

because the curves shown in Fig.Q]are almost identical up 
to an angle scaling. We cut the distribution at 6 = tt/2 
to avoid backward-going modes. We assume monoener- 
getic neutrinos with lo c = Am 2 /2E C — 0.5 km -1 for the 
core flux, i.e. E c ~ 12 MeV. We have checked that this 
choice reproduces well the instability region in Fig. 4 of 
our previous paper [8J where the full spectrum was used. 
The halo energies are larger as described in Sec. Ill Bl im- 
plying (E^ 1 ) = 2(E^ 1 ), thus we use monoenergetic halo 
neutrinos with u)h = 0.25 km - . 

The system is unstable in flavor space if the linearized 
equation of motion discussed in Ref. |8J has eigenvalues 
with imaginary part k, the radial growth rate (in units of 
km^ 1 ) of the unstable mode. Whether this growth rate is 
large or small depends on the available distance. In other 
words, at a given distance r, the instability is important 
if nr 3> 1 and unimportant if nr <C 1. Therefore, in Fig. [5] 
we show contours of constant nr in a two-dimensional pa- 
rameter space consisting of radius r of our model and an 
assumed electron density n e which causes the multiangle 
matter effect. 

If we use the core flux alone, this figure contains the 
same information as Fig. 4 in our previous paper Q. In 
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FIG. 6: Contours for the indicated nr values in the parameter 
space of radius r and assumed electron density n e . We use 
a simplified energy spectrum as described in the text. The 
core flux alone is responsible for Kr > 1 for r < 700 km (blue 
shading), the halo adds the red-shaded region at larger r. We 
also show the density profile of our 280 ms numerical model. 

the absence of matter, the instability would set in at 
about 150 km, whereas the presence of matter suppresses 
the instability entirely: the SN density profile never over- 
laps with the instability region. The most dangerous re- 
gion of closest approach is at around 600 km. 

In Fig. |6]one clearly sees how the halo flux extends the 
instability region to larger radii (red shading), but only 
far away from the SN density profile. In other words, the 
multiangle matter effect strongly suppresses the would- 
be instability caused by the halo. 

This finding does not contradict the importance of the 
halo flux for neutrino-neutrino refraction. It is true that 
the neutrino-neutrino interaction energy caused by the 
halo was found to decrease roughly as r -3 , instead of r~ 4 
by the core flux alone. However, the multiangle matter 
effect is enhanced in a similar way because it depends on 
the same neutrino angular distribution. A more appro- 
priate comparison is between the neutrino and electron 
densities. The former decreases as r~ 2 whereas in our 
model the latter decreases roughly as r -135 and thus be- 
comes relatively more important with distance. There- 
fore, it is no surprise that the multiangle matter effect 
would be more important with increasing distance. 



IV. CONCLUSIONS 

The halo distribution of neutrinos causes neutrino- 
neutrino refraction to decrease more slowly with distance 
than had been thought, but still faster than the matter 
density during the accretion phase. Repeating our lin- 
earized stability analysis including the halo, the system 
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remains stable for the chosen example. As anticipated, 
the multiangle matter effect is a crucial ingredient for 
collective flavor conversions especially in those regions 
where the halo flux is important. 

Thus it remains possible that self-induced flavor con- 
versions are generically suppressed during the early ac- 
cretion phase. The early signal rise would then encode a 
measurable imprint of the neutrino mass hierarchy fTlj | . 

However, more systematic studies are required to un- 
derstand if such conclusions are generic. For example, 
we have assumed that the relatively small backward flux 
can be neglected for the stability analysis, but its possible 
role needs to be better understood. Likewise, the role of 
multi-D effects in SN modeling remains to be explored. 

The topic of collective flavor oscillations of SN neutri- 
nos remains a fiendishly complicated problem. 
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Appendix: Analytic Halo Estimate 

In order to derive the analytic halo estimate of Eq. (|4| 
we consider a radial distribution of scattering targets, 
n(£), and a differential scattering cross section da/dQ. 
The scattering rate per volume and per solid angle for 
neutrinos with number density n v {£) at radius £ is 



dq s , 



dtt 



n v {£)n(£)c 



da 

~dn' 



(A.l) 



where we have reinstated the speed of light. The "number 
intensity" I(r, 9) at radius r of neutrinos scattered to 
angle direction 9 relative to the radius vector at r can be 
obtained by adding up all scattering contributions along 
the path s as shown in Fig. [7] It depends on the total 
neutrino emission rate $ = Lj (E) of the pointlike source 
assumed to be located at r — £ = 0, which determines 
the neutrino number density n u (£) = $ /(An£ 2 c). With 
this result and Eq. (| A. 1 1) one finds 



X(r, 



dn 



da 



ds 4^ n( - e) dn ( - a + l 



(A.2) 



FIG. 7: Geometric definitions for the analytic halo estimate. 



where d = a + 8 is the scattering angle and £ = 
r sin 8 1 sin $ and s = r sin a/ sin $, see Fig. [7J 

Using a decreasing power-law distribution of the target 
density, n(£) — nn{R/£) m with m > 0, an assumed de- 
pendence of the scattering cross section on the scattering 
angle § of the form da($)/dfl = <7 S catt(l + acos$)/47r, 
one can transform the integral of Eq. (|A.2[) into the form 
of Eq. (4). 
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